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SUMMARY 


A three-dimensional subsonic aerodynamic panel code (VSAERO) was used to 
predict the effects of upper and lower external nozzle flap geometry on the external 
afterbody/nozzle pressure coefficient distributions and external nozzle drag of 
nonaxi symmetric convergent-divergent exhaust nozzles having parallel external 
sidewalls installed on a generic twin-engine high performance aircraft model. 

Nozzle static pressure coefficient distributions along the upper and lower surfaces 
near the model centerline and near the outer edges (corner) of the two surfaces were 
calculated, and nozzle drag was predicted using these surface pressure 
distributions. A comparison between the theoretical predictions and experimental 
wind tunnel data is made to evaluate the utility of the code in calculating the flow 
about these types of non-axi symmetric afterbody configurations. 

For free-stream Mach numbers of 0.60 and 0.90, the conditions where the flows 
were attached on the boattails yielded the best comparison between the theoretical 
predictions and the experimental data. For the boattail terminal angles of greater 
than 15°, the experimental data for M = 0.60 and 0.90 indicated areas of separated 
flow, so the theoretical predictions failed to match the experimental data. Even 
though calculations of regions of separated flows are within the capabilities of the 
theoretical method, acceptable solutions were not obtained. 


INTRODUCTION 

Recent studies (ref. 1-12) of nonaxi symmetric nozzles have shown potential 
improvements In aircraft performance and effectiveness especially with the 
incorporation of thrust vectoring and reversing into the nozzle concept. It is also 
understood from these references that once non-axi symmetric nozzles are incorporated 
into an aircraft concept, thrust vectoring and thrust reversing can be added with 
only a minor weight penalty. Thrust vectoring and reversing, by virtue of its 
potential for enhanced aircraft capability, may be a requirement for the next 
generation fighter aircraft. If these nozzles are to be utilized they must be 
integrated into the aircraft in such a way as to either reduce the aircraft drag or 
maintain the drag at levels associated with current axi symmetric nozzles. A number 
of studies have been undertaken to determine the most efficient method for 
integrating non-axi symmetric nozzles into twin engine fighter aircraft. The 
development and utilization of advanced computational methods will play a vital role 
in developing this technology. As a result the current form of research activities 
at the Langley Research Center is to develop the data base for new non-axi syimetric 
nozzle concepts, as well as, to develop and validate computational methods for 
predicting the flow characteristics on these nozzles. 

The results of an investigation to determine the external drag characteristics 
of non-axi symmetric twin engine type afterbody /nozzle shape has been reported in 
reference 13. In this reference the effects of upper and lower external nozzle flap 
geometry on the afterbody and external nozzle pressure coefficient distributions, 
afterbody forces and moments, and integrated nozzle axial force are presented. The 
nozzle flap geometry was formed by a parametric combination of circular arcs and 
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straight lines with the nozzle external sidewalls having parallel surfaces. For 
this report, the experimental data for configurations 2, 5, 13, 14, and 21 at M = 
0.60 and M = 0.90 at an angle-of -attack of zero degrees (a = 0°) and operating 
nozzle pressure ratios (NPR) from jet-off to 8 are compared to the theoretical 
predictions of the vortex separation aerodynamic program (VSAER0) detailed in 
references 14 and 15. These configurations offer a comparison between nozzles of 
the same length with varying terminal boattail angle, nozzles of the same terminal 
boattail angle with varying circular arc fairing radii, and one very short steep 
boattail profile. 

VSAER0 is a surface singularity program using quadrilateral panels to represent 
arbitrary three-dimensional bodies. This code includes such features as wake shape 
iteration, jet-wakes, on and off body stream line calculations, viscous/inviscid 
interactions and off-body velocity calculations. Other features include modeling of 
jet exhaust and separated wakes and calculating the effects of flow separation and 
boundary layer growth. This code was used to calculate the external static pressure 
coefficient distributions and nozzle drag coefficients for the subject 
nozzle/afterbodies. The results are then compared to the experimental data of 
reference 13 to assess the capability of the code to predict 2-D C-D 
afterbody /nozzle external flows and external drag. 


SYMBOLS 


Model forces and moments are referenced to the stability-axis system with the 
model reference center located at 0.25c, which corresponds to fuselage station 
{ FS ) 36.09 Inches and 1.75 inches above the model centerline. All aerodynamic 
coefficients are ndndlmenslonalized with respect to q^S or q M Sc . 

Cq drag coefficient, D/q^S 

C D>aft afterbody drag coefficient 

c D,n 


c 


nozzle drag coefficient 

P-P. 

pressure coefficient, — 

Am 

wing mean geometric chord, 17.47 in. 


M, MACH 
P 

P 

00 

R 

S 

W 


length, in. 

Mach number 

surface static pressure, psl 
free-stream static pressure, psi 
free-stream dynamic pressure, psi 
radius, in. 

p 

wing reference area, 664 in c 
width, in. 
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x, X 

y 

a, ALPHA 

e 

Subscripts: 

A 


c 

e 

f 

1 

n, N 

o 

t 

Abbreviations 

A/B 

BL 

Conf . 

FS 

NPR 

WL 


distance downstream from FS 66.54 
spanwise distance from model centerline, in. 
model angle of attack, deg 
boattail angle, deg 

afterbody 

chord 

nozzle exit 
flap 

arc/straight line tangent point 
nozzle 

circular arc center 

throat location or terminal angle 

afterburning 

butt line 

configuration 

fuselage station 

nozzle pressure ratio, p t j/p^ 

water line 

two-di mensi onal convergent-di vergent 


EXPERIMENTAL DATA AND THEORETICAL METHOD 
Wind Tunnel 

The experimental investigation of reference 13 was conducted in the Langley 16- 
Foot Transonic Tunnel, a single-return, continuous flow, atmospheric-pressure wind 
tunnel with a slotted, octagonal test section and continuous air exchange for 
cooling. The wind tunnel has a variable airspeed up to a Mach number of 1.30. Test 
section plenum suction is used for speeds above a Mach number of 1.05. A complete 
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description of this facility, its related equipment, and its operating 
characteristics can be found in references 16 and 17. 


Model and Tests 

A sketch of the wing-tip supported model with a typical nozzle installed is 
shown in figure 1(a) with details of the jet simulation system and balance 
arrangement. The overall model arrangement represents a twin engine fighter with 
tail surfaces removed and faired over, twin side inlets. Figure 1(b) is a half span 
sketch showing the wing planform geometry and model reference moment center. The 
wing was mounted 1.75 inches above the model centerline, had an aspect ratio of 
2.40, a taper ratio of 0.43, and a NACA 64-series airfoil section with thickness 
ratio of .067 near the root. The outboard area of the wing had increased thickness 
ratio (0.10 at the tip) for structural support and for compressed air and 
instrumentation passages. The wing/centerbody together with the twin support booms, 
formed the bifUrcated-stlng model support system. A complete description of this 
support system and the various balance arrangements that may be used can be found in 
references 6 and 16, and further description of the arrangement used in these 
particular tests may be found in reference 13. 

Photographs of the model installed in the test section of the Langley 16-Foot 
Transonic Tunnel are presented in figure 2. 

The twin-eriglne afterbody shown in figure 3(a) and 3(b) was designed to be 
representative of current fighter aircraft, while housing the propulsion simulation 
system, afterbody balance, and related static pressure tubes from the nozzles and 
afterbody. The afterbody is constant in width so that all closure is accounted for 
on the top and bottom surfaces using a 5° boattail between FS 59.19 and the end of 
the afterbody at FS 66.04 (see figure 3(a)). 

Figure 4(a) presents a sketch and table giving geometry of the two-dimensional 
convergent-divergent ( 2— D C-D) nozzle models tested. Of the 21 configurations shown 
in the table arid repdrted in reference 13, only the 5 indicated by the + sign in the 
configuration coliimh are used in this paper. The five configurations selected are 
dry power corifigUrations using identical internal nozzle geometry. 

An external High-pressure air supply system provided the nozzles with a 
continuous flow of clean dry air at a controlled temperature of 70°F. The air was 
used to simulate exhaust flow over a range of nozzle pressure ratios (NPR's) from 
about 1.0 (jet off) Up to about 10.0 depending on Mach number. 

Figure 4(b) shows the locations for the 23 static pressure orifices on the top 
flap, and the 11 orifices on the bottom flap of each nozzle. Each orifice number 
was located at the same ratio of nozzle length (x/l n ), regardless of nozzle 
configuration. 

The external static pressure orifice data on the right-hand nozzle were used to 
compute external forces and moments on the entire nozzle by a numerical summation of 
the local nozzle static-pressure coefficient multiplied by a projected area assigned 
to each of the orifices on the nozzle and all divided by the reference area. Flow 
on the model was assumed to be symmetrical right to left, so each orifice was 
assigned both right and left side areas appropriate to the particular orifice. 

Normal and axial projected areas were used to compute drag coefficient data. 
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Model angle of attack was computed from an attitude indicator located inside of 
one of the boom fairings, and corrected for 0.1° flow angularity, which is the 
average flow angle measured in the 16-Foot Transonic tunnel. 

The investigation of reference 13 was conducted in the Langley 16-Foot 
Transonic Tunnel at free-stream Mach numbers of 0.60, 0.70, 0.80, 0.90, and 1.20. 

The nozzle pressure ratio was varied from 1.0 (jet off) to 10.0, depending on free- 
stream Mach number. For this report, only data for M = 0.60 and M = 0.90 at 
a = 0° are used with NPR set at typical operating values for each Mach number. 


VSAERO Theoretical Prediction Method 

VSAERO (Vortex Separation Aerodynamic Program) is a surface singularity 
solution to the Laplace equation using quadrilateral panels to represent arbitrary 
three-dimensional bodies. Source and doublet singularities are distributed in a 
piecewise constant fashion on each panel. Neumann boundary conditions are applied 
to the panel to determine the singularity strength. Dirichlet boundary conditions 
are applied at the boundary to determine the doublet strength. More detailed 
discussions of the method appear in references 14 and 15. 

The code includes features such as wake shape iteration, jet exhaust, on- and 
off-body streamline calculations, and off-body velocity calculations. Jet wakes, in 
particular, are modeled using doublet sheets with a linear variation in doublet 
strength in the streamwise direction. Jet exhaust flow can be modeled by specifying 
a unit normal velocity outflow at the exhaust. The major limitation in the use of 
the code, is the incompressible nature of the basic equation. However, transonic 
results may be approximated through the use of the Prandtl-Glauert compressibility 
correction. The VSAERO code is capable of modeling the full geometry of the 2-D C-D 
model configurations including body, wing, exhaust nozzle, and wing-tip support 
booms. 

The geometrical description of the various configurations is based on panels, 
patches, and wakes, which can be combined to form components and assemblies. Some 
automatic capabilities are provided in the code for paneling and also patch 
generation in special situations. The computed geometrical description of the 2-D 
C-D afterbody model is shown in figure 5 as produced from data input to the VSAERO 
code. 


The effects of viscosity are calculated using a two-dimensional, integral 
boundary layer calculation (ref. 18). The calculations follow surface streamlines 
that have been determined from an initial potential flow calculation. The boundary 
layer thickness from this calculation is simulated by surface transpiration, and a 
new potential calculation is performed using the transpiration velocities. 

Typical computational times for the 2-D C-D boattail model using approximately 
1100 panels and running 6 to 8 boundary layer iterations were 4 hours on a MicroVAX 
II. 
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RESULTS AND DISCUSSION 


Nozzle Pressure Distribution 


Static-pressure-coefficient distributions on the external surface of the nozzle 
boattail are presented in figure 6 for M = 0.60 data and figure 7 for M = 0.90 data 
for configurations 2, 5, 13, 14, and 21 of table 4(a). All data are for a zero 
angle-of-attack and symbols are used to indicate the experimentally determined 
nozzle-pressure distributions. The location of the rows of orifices on the boattail 
are shown in figure 4(b). Theoretical pressure distributions computed are presented 
in figures 6 and 7 as solid lines. 

Comparisons at M =0.60 .- In general, the agreement of the theoretical and 
experimental pressure distributions are good (see figure 6). Typical nozzle 
boattail pressure distributions are characterized at subsonic speeds by an expansion 
near the start of the boattail followed by a pressure recovery on the remaining 
portion of the boattail (x /I = 0) • The degree of flow expansion depends upon 
the rate of change cif slope n between the model afterbody and the nozzle boattail at 
the nozzle connect station. 


The degree of pressure recovery depended Upon the nozzle pressure ratio, free- 
stream Mach number and, to some degree, upon proximity to the nozzle external 
sidewall. The theoretical method predicted the degree of initial flow expansion for 
configurations in which the slope transition from afterbody to nozzle boattail was 
small or a circular-arc fairing (see figure 6(a) for example). Little or no flow 
separation at the start of the boattail, i.e., Xg/l n = -0.2 to +0.2, was evident for 
eveh the shortest ahd steepest nozzle, configuration 5. Poor agreement occurred at 
the nozzle cbrrter hear the sidewall, i.e., rows 4 and 5 for all nozzle concepts. 
Also, poor agreement occurred at the nozzle trailing edge where jet entrainment and 
plume blockage effects tend to dominate the flow. For these calculations, the shape 
of the simulated jet Wake was allowed to relax to assume a shape determined by the 
pressures calculated by the initial potential flow solutions. This new shape was 
used for any subseqlieht viscous and inviscid flow solutions. The jet wake 
relaxation colild account for some part of a plUme blockage effect. Jet entrainment 
effects were not included in these calculations, though a type of entrainment 
modeling Was Within the capability of the code. The predicted spanwise variation of 
the level of pressure recovery, by comparing the point x n /lp = .9 for rows 1 through 
4, was greater than that for the experimental data. This slight mismatch in the 
comparison is caused by a lack of accurate plume modeling. The poorest agreement, 
overall, occurred for configuration 5 (see figure 6(b)) in the region of apparent 
flow separation, x n /1p = 0.5 to 1.0 for rows 1 and 2. Though separated wake 
modeling was also a capability of the code, a reasonable and stable solution 
utilizing that modeling was not obtained. This modeling problem was likely due to 
insufficient damping of the flow solution in the separated flow region between each 
boundary-layer iteration. In addition, tile area of flow separation appeared 
irregular in shape, i.e., not along a fixed row or column of the paneled geometry, 
thus subsequently making a good separated flow solution difficult to obtain. 


Changes in the pressure coefficient distribution with variation in nozzle 
cross-section shape were predicted well. Abrupt changes in surface curvature were 
particularly noticeable in the configurations with combination circular arc-straight 
line shape cross-sections. Configuration 14, figure 6(d) and 7(d), at x n /l n = .4 
and configuration 21, figure 6(e) and 7(e), at x n /l n = 0.05, both show 
discontinuities in the pressure coefficient distributions due to change in surface 
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curvature. Though it was difficult to assess the accuracy of the predicted 
magnitude of the expansion peak, the experimental data bracketing this region are 
matched fairly well. 

Comparisons at M = 0.90 .- Overall, the observations made for the comparisons 
at M = 0.60, held for M = 0.90 (see figure 7) due to the external flow expansion, 
and following recompression remained subsonic and attached for all configurations 
except configuration 5. All solutions computed for M = 0.90 were for inviscid flow 
due to convergence problems with the boundary-layer iteration scheme at this Mach 
number. The slope discontinuity at the connect station of configuration 5 (fig. 
7(b)) was sufficient to cause flow separation off the surface of the boattail. 
Therefore, there was poor agreement of the predicted pressure coefficients with the 
experimental data. The inviscid attached flow solution greatly over-predicted the 
initial flow expansion about the start of the boattail and over-predicted the 
subsequent flow recompression. The largest discrepancy between the theoretical and 
experimental data occur in the region of the nozzle trailing edge where jet and 
viscous effects are predominant. 


Drag Comparison 

Theoretically predicted afterbody /nozzle drag coefficient, and experimental 
drag coefficient for M = 0.6 and M = 0.90 are compared in figure 8. The discussion 
will be limited to that of the 10° nozzle terminal boattail angle configurations; 2, 
14, and 21. The major difference in the boattail longitudinal shape between these 
three configurations was the radius of curvature of the circular arc, fairing the 
flat portion of the nozzle to the afterbody, i.e. column heading R in table 4(a). 

As indicated by the sketch at the top of figure 8, nozzle drag coefficient is 
indicated by an open bar and afterbody drag coefficient by a cross hatched bar, and 
the total afterbody /nozzle drag coefficient by the arrow. Configuration number is 
indicated at the bottom of each bar. As you can see from these data the integrated 
pressure distributions on the afterbody and nozzle do not agree with the measured 
drag data. This can be attributed to the fact that the pressure drag was calculated 
by the area-pressure integration procedure which always yields some inaccuracies 
simply because an adequate number of measurements is difficult to obtain. 

The prediction of the lowest afterbody drag, configuration 14, agreed with the 
experimental test data for both Mach numbers, though the absolute level of drag was 
considerably underpredicted. The predicted trend of total afterbody /nozzle drag 
coefficient for the three configurations did not agree with the test data at either 
Mach number. Different trends in the nozzle drag coefficient occur at M = 0.6 
comparing the theory with experimental data. At M = 0.9, the degree of change in 
afterbody drag coefficient due to changing nozzle shape was sufficient to cause the 
disagreement in the total drag trends. A partial reason for this discrepancy can be 
explained by the fact that the experimentally determined nozzle pressure drag was 
calculated by the area-pressure integration procedure. The arrangement of the 
orifices used for the integration and the number of orifices used can introduce some 
error into the calculation. For this set of nozzles, there were no pressure data 
obtained for the outboard 10 percent of the nozzle planform. It is interesting to 
note that the theoretical predictions show the corner region to be an area of 
positive drag, as shown in figure 9. This figure is a plot of integrated pressure 
drag along lines of constant span, moving from the nozzle boattail centerline to the 
sidewall. It appears that the central 90 percent of the boattail surface 
contributes a negative drag component to the total drag of the nozzle. The strong 
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corner flow turning from the sidewall onto the upper and lower external nozzle flaps 
induce negative pressures in this region, resulting in the positive drag 
component. The magnitude of this additional drag did not change significantly 
between the three configurations, and may account for the differences In the level 
of drag of each nozzle (between theory and the experimental data) but not the trend 
of drag change between each nozzle. 


CONCLUSIONS 

A three-dimensional subsonic aerodynamic panel code (VSAERO) was used to 
predict the effects of external nozzle flap geometry on the external 
afterbody /nozzle drag of nonaxi symmetric convergent-divergent exhaust nozzles having 
parallel external sidewalls installed on a generic twin-engine high performance 
aircraft model. External geometries for 5 configurations were used to evaluate the 
utility of the code. A comparison of the static pressure coefficients predicted on 
the nozzle boattail with the experimental data and a comparison of afterbody /nozzle 
drag coefficients result in the following conclusions: 

1. Good prediction of the flow about afterbody/nozzle configurations 
can be made, except in cases of separated or supersonic flows. The 
computational code did predict the occurrance of flow separation, however 
a converged solution of the separated flow was not obtained. 

2. Prediction of total afterbody /nozzle drag trends did not agree with 
the experimental data. The incremental changes in drag coefficient 
for both the experimental data and the theoretical predictions were 
small. A clear advantage of one configuration over another was difficult 
to discern. 
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Section A-A 


(a) Jet simulation system and balance arrangement. 
Figure 1.- Sketch of air powered, twin engine, wing-tip supported 
model with two-dimensional convergent-divergent 
nozzles from reference 13. All linear 
dimensions in inches. 







DOWNSTREAM VIEW 


Figure 2.- Photographs fr 
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Figure 3.- Afterbody geometric description and static pressure orifice 
locations from reference 13. All linear dimensions 
-■ in inches 
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Figure 3.- Concluded 
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♦ Configurations used in this comparison. 

(a) Nozzle geometry. 

Figure 4.- Sketches of nozzle assembly, geometric design parameters, and 
static pressure orifice locations on the nozzle, from 
reference 13. All linear dimensions in inches. . 
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Figure 5.- 


1 mnnf? of paneling of wing-tip supported afterbody 
model used for input to the Vortex Separation y 
Aerodynamic Program (VSAERO). 







(a) Configuration 2. 

Figure 6.- Surface static pressure coefficient distributions around the nozzles 
at M = 0.60 and a = 0°. Solid lines represent the computed values. 
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(c) Configuration 13. 


Figure 6.- Continued 
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(a) Configuration 2. 

Figure 7.- Surface static pressure coefficient distributions around the nozzles 
at M = 0.90 and a = 0 . Solid lines represent the computed values. 
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CONFIGURATION 5. 
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(b) Configurat 
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(c) Configuration 13. 


Figure 7.- Continued 
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Figure 8.- Comparison of theoretically predicted afterbody/nozzle 
drag coefficient and experimental drag coefficient for 
8 t = 10°, at M = 0.60 and 0.90, and a = 0°. 
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